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ABSTRACT 

A Fourier analysis on galaxy number counts from redshift data of both the Sloan Digital 
Sky Survey and the 2dF Galaxy Redshift Survey indicates that galaxies have preferred periodic 
redshift spacings of Az = 0.0102, 0.0246, and 0.0448 in the SDSS and strong agreement with 
the results from the 2dF GRS. The redshift spacings are confirmed by the mass density fluctua- 
tions, the power spectrum P{z) and Npairs calculations. Application of the Hubble law results in 
galaxies preferentially located on co-moving concentric shells with periodic spacings. The com- 
bined results from both surveys indicate regular co-moving radial distance spacings of 31.7 ± 1.8 
h~^Mpc, 73.4 ± 5.8 h~^Mpc and 127 ± 21 h~^Mpc. The results are consistent with oscillations 
in the expansion rate of the universe over past epochs. 

Subject headings: galaxies:distances and redshifts — large-scale structure of universe 



1. Introduction 

When modeling the large scale structure of 
the cosmos the cosmological principle is assumed, 
therefore what we see must be biased by our view- 
point. The assumption is that the universe has 
expanded over time and any observer at any place 
at the same epoch would see essentially the same 
picture of the large scale distribution of galaxies 
in the universe. 

However evidence has eme rged from both the 
2dF Galaxy Reds hift Su rvey (l2dF GRSl l (see fig. 
17 of CoUess et al ('200l|)) and also the Sloan Dig- 
ital Sky Survey (SDSS) (see fig. 2 of Tegmark et 
al (|2004D 1 that seem to indicate that there is peri- 
odicity in the abundances of measured galaxy red- 
shifts. This has emerged from the galaxy number 
counts (iV) within a small redshift interval {6z) as 
a function of redshift (z). 

What has been found is the so-called picket- 



fence structure of the N-z relation, f irst n oticed 
in one dimension by Broadhurst et al ( 1990l ) from 
a pencil-beam survey of field galaxies. In this pa- 
per the galaxy surveys are analysed on 2D and 
3D scales. The usual interpretation is that this 
is evidence of the Cosmic Web, including voids 
and long filaments of galaxies. An alternative in- 
terpretation is suggested in this paper; galaxies 
are preferentially found in redshift space at cer- 
tain redshifts with higher number densities than 
at other redshifts. Of course, this does not rule 
out cosmic web structures in addition to this ef- 
fect. 

However the concept that this represents real 
space galaxy distances is incompatible with the 
cosmological principle, which assumes the unifor- 
mity of space at all epochs on sufficiently large 
enough scales. And it has been demonstrated 



([Yoshida et al 20011 ). from a large number of nu- 
merical simulations using the Einstein-de Sitter 
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and ACDM models, that the probabiHty of get- 
ting such a periodic spatial structure from cluster- 
ing and Cosmic Web filaments is less than 10~^. 
This allows for an alternate explanation, namely 
that the redshift periodicity is evidence for past 
oscillations in the expansion rate of the universe. 

The only information that has been sought from 
large scale galaxy surveys is the spatial power 
spectrum that describes the random yet uniform 
distribution of galaxies in the universe, where the 
only departure expected from the random distri- 
bution is some characteristic scale for the cluster- 
ing of galaxies. Certainly very few expect to find 
evidence consistent with a periodic redshift distri- 
bution. 

However Hirano et al (l2008h have attempted to 
explain galaxy redshift abundance periodicity as 
a real effect resulting from the universe having a 
non-minimally coupled scalar field that manifests, 
in later cosmological times. The effect produced 
oscillations in the expansion rate as a function of 
time. They model not only the redshift space ob- 
servations but also the effect their cosmology has 
on the CMB power spectrum and the type la su- 
pernova Hubble diagram. 

In the usual analysis the spatial two-point or 
autocorrelation function is used to define the ex- 
cess probability, compared to that expected for a 
random distribution, of fin ding a pair of galax- 
ies at a given separation ( Baugh 20061 ). As a 
result the power spectrum P{k) i s derived from 
the two-point correlation function (jPeebles 
The power spectrum is predicted by theories for 
the formation of large scale structure in the uni- 
verse and compared with that measured, or more 
precisely calculated from the available data. 

The power spectrum for the SDSS has been 
calc ulated using a set of basis functions defined 
in llHamilto n et al. 2000l ). Three power spec- 
tra (galaxy-galaxy, galaxy-velocity, and velocity- 
velocity) were determined. Local motions of 
galaxies only contribute a radial component of a 
galaxy's total redshift, hence only affect the radial 
component of the power spectrum. Because the 
angular power spectrum is unaffected the assump- 
tion is also made that galaxy-galaxy spectrum is 
equal to the redshift-space power spectrum in the 
transverse direction. In this paper we analyze 
the redshift hence radial component of the power 
spectrum P{z), calculated in redshift space not 



fc-space, and find periodicity consistent with that 
observed in N(z). Also an Npairs analysis is car- 
ried out on both survey data sets with confirming 
results. 



The "Bull's-eye" e ffect has been studie d (|Praton et al. 1997 
iMelott et al. 1998t iThomas et al 2004 via A^- 
body simulations to result from large-scale in- 
fall plus small-scale virial motion of galaxies. It 
is believed that these two effects can bias mea- 
surements of large scale real space galaxy dis- 
tributions. One is the so-called "fingers-of-God" 
(FOG) effect, which results from random peculiar 
motions of galaxies at the same distance from the 
observer. The second acts on much larger scales 
and where overdensities of galaxies occur, like at 
the "Great Wall" for instance. See Fig. 1. 

Galaxies tend to have local motions toward 
the center of such structures. These motions are 
not random but coherent and add or subtract 
to the observer's line-of-sight redshift determina- 
tion. These effects preferentially distort the map 
in redshift space toward the observer due to the 
velocities of galaxies within clusters, i.e. non- 
cosmological redshift contributions. However the 
latter can only enhance in redshift space existing 
weak real space structures, it cannot create con- 
centric structures centered on the Galaxy. And 
the FOG effect tends to smooth out finer detailed 
or higher Fourier frequency structures in redshift 
space and is not very relevant to this discussion. 

In the maps, fig. 18 of Colles s et al (|200ll) ) 
and fig. 4 of Tegmark et al (|2004l ). there appears 
to the eye to be concentric structure in the ar- 
rangement of galaxies. See the left hand side of 
the Fig. 1, the SDSS data with declinations be- 
tween 0° and 6° mapped in Cartesian coordinates. 
Could it be that we inhabit a region of the Uni- 
verse that is a bubble, part of a background of 
many such bubbles? But because redshift only 
gives us one dimensional information the conclu- 
sion could be drawn that if the galaxies in the 
Universe are evenly distributed throughout then 
redshift periodicity is indicative of oscillations in 
the expansion at past epochs. All observers would 
see the same shell structure in redshift space re- 
gardless of their location. 

In this analysis we take a heuristic approach 
and look for periodicity in the redshift spacing be- 
tween galaxies, in 2D in the case of slices and in 
3D in case of the whole SDSS data set. In practice 
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the discrete Fourier Transform {vs) is calculated 
from the N-z relations, which are histograms de- 
termined by binning (counting) the observed red- 
shifts of the survey galaxies between z — 6z/2 and 
z + 5z/2, as a function of redshift z. Within the 
ith bin there are ANi galaxies. Discrete Fourier 
amplitudes Vs are generated for each integer s > 1 
from 

:.. = ^Vuse2'^^"(-i)(^"i)/", (1) 

where Ug = ANi, j = 1, is the total number 
of ANi bins. Because are complex their abso- 
lute value is taken in the analysis. The Fourier 
frequencies are calculated from 



where Az is the periodic redshift interval. 

The results suggest that there are redshifts, on 
the scale size discussed here, which is not primar- 
ily the result of clustering or the FOG and infall 
effects resulting in distortion of the redshift space 
maps. Both the SDSS and 2dF GRS redshift data 
where the redshifts are known with 95% confidence 
or greater are used. Fourier analysis can find pe- 
riodic structure even where it may not be so ap- 
parent to the eye. 

Data for 427,513 galaxies from the SDSS Fifth 
Data Release (DR5) were obtained where the data 
are primarily sampled from within about -10 to 
70 degrees of the celestial equator. Also data for 
229,193 galaxies were obtained from the 2dFGRS 
where the data are confined to within 2 degrees 
near the celestial equator and balanced between 
the Northern and Southern hemispheres. In the 
SDSS case the data are not so well balanced, yet 
there are nearly twice as many usable data. 

2. Fourier analysis of N-z relations 
2.1. 2D slices of SDSS data 

From the SDSS data three 6° slices were taken 
in declination (Dec) angle for examination and 
comparison. They were taken from 0° < Dec < 6° 
, 40° < Dec < 46° and 52° < Dec < 58°. A real 
space Cartesian plot (with aspect ratio of unity) 
of the data from the first slice near the equato- 
rial plane is shown in Fig. 1. The redshifts have 



been converted to co-moving radial distance with 
Hubble constant Hq = 100 kms~^Mpc~^ where 
the parameter h = i/o/100. The Great Wall is 
shown in the second and third quadrants as indi- 
cated. In those two quadrants it is evident to the 
eye that there is general concentric structure with 
a spacing of about 75 h^^ Mpc. 

Polar coordinate maps of these regions, pro- 
jected onto the equatorial plane, are shown in Figs 
2, 3 and 4, respectively. Most of the SDSS data 
come from two declination bands on the sky. The 
first slice is taken from the lower band near the 
celestial equator and the latter two from higher 
declinations. Polar maps can tend to distort a real 
space image but a comparison of Figs 1 and 2 in- 
dicates that these polar plots would be very close 
in appearance to real space maps if converted to 
co-moving coordinates. 

From simple inspection concentric structure is 
apparent in both Figs 2 and 4 but not as apparent 
in Fig. 3. To the eye this seems too coinciden- 
tal and deserves further investigation. At least at 
some declinations the case can be made for red- 
shift periodicity. 

These data were binned with a bin size 5z — 
10~^ resulting in a N-z relation. The redshift res- 
olution of most, but not all, of the SDSS data 
is only 10~^. This bin size essentially counts all 
galaxies of the same redshift regardless of their 
position on the sky. 

Using ll]) the discrete Fourier spectra were cal- 
culated for N{z) of each data slice, and are shown 
in Fig. 5. The Fourier frequency axis in Fig. 5, 
and all subsequent figures, has been converted to 
a redshift interval (Az). For the lowest angle slice 
(solid black curve) from the data of Fig. 2 and 
the highest angle slice (solid gray curve) from the 
data of Fig. 4 the peaks at numbers 1 and 2 co- 
incide. For the data from the 40° < Dec < 46° 
slice of Fig. 3 the peaks are not strong nor do 
they coincide with the others except at peak num- 
ber 1. The fact that the maps in Figs 1 through 3 
are visually different and that the resulting Fourier 
spectra don't all have coincidental peaks indicates 
that the peaks that do coincide are not the result 
of some selection effect resulting from the choice 
of filters or otherwise. Therefore it is evident that 
the peaks are the result of real redshift structure 
that varies in different directions. The underlying 
cause must be common to all, since all three slices 
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Fig. 1. — Real space Cartesian coordinate plot 
in units of Mpc: co- moving (Hubble) distance 
derived from SDSS redshift data from 0° < Dec < 
6°. 
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Fig. 3. — Polar coordinate plot as a function of 
RA: Redshift space map of data from 40° < Dec < 
46°. There are 44,335 data for z < 0.5. There are 
no observed data at these declinations in the two 
empty quadrants. 



90 








240 300 
270 



Fig. 2. — Polar coordinate plot as a function of 
RA: Redshift space map of the same data as Fig. 
1, taken from 0° < Dec < 6°. There are 59,930 
data for z < 0.5. 
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Fig. 4. — Polar coordinate plot as a function of 
RA: Redshift space map of data from 52° < Dec < 
58°. There are 34,802 data for z < 0.5. There are 
no observed data at these declinations in the two 
empty quadrants. 
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have a Fourier peak at number 1, which represents 
a redshift space interval of Az w 10~^. 

2.2. 3D: All SDSS data 

In order to see if large scale 3D concentric struc- 
ture is evident the whole SDSS galaxy data set was 
binned and Fig. 6 shows the resulting N-z relation 
for all z < 0.35. The discrete Fourier spectrum 
was calculated from N{z) using ([T]) and the result 
is shown in Fig. 7. It was found that z > 0.15 
have very little effect on the result. 

Also the SDSS data were divided into two sets 
sampled from two opposite regions of the sky. The 
data are 374,425 with right ascensions 100° < 
RA < 270°, and 53,087 with 290° < RA < 50°. 
There is about 7 times the number of data in one 
set than the other. The larger set is mostly from 
above the celestial equator. 

N(z) relations for these regions are shown in 
Fig. 8 on a double-Y axis since they are about 
7 times different in scale. Nevertheless the gen- 
eral shape of the curves are the same indicating 
that N{z) has the same z-dependence in each re- 
gion. The location of peaks on these curves are 
similar except peaks occur with different relative 
heights and some at different redshifts. This indi- 
cates that though N{z) has the same general de- 
pendence (shape of the galaxy survey function) the 
different regions are distinctly different in other 
respects. Again this is evidence that the N-z rela- 
tions are not artifacts of the filters or due to some 
random effects. 

This is seen also in their respective Fourier 
transforms as shown in Fig. 9. Curve 1 is the 
result from the whole SDSS data set (as shown 
in Fig. 7), curve 2 is the Fourier spectrum from 
the larger region of the sky surveyed with 100° < 
RA < 270° and curve 3 (enhanced by a factor of 
3 for better comparison) is from the smaller re- 
gion of the sky with 290° < RA < 50°. There are 
similar features but much weaker in the latter, yet 
also it is apparent that the peaks are shifted in 
Fourier frequency. This could indicate real differ- 
ences in the redshift space structure in that part 
of the survey as compared to the other. 

In order to evaluate the level of significance of 
the peaks from the whole SDSS data set the follow- 
ing was done. To the galaxy survey function, the 
smooth polynomial that fits the N-z data of Fig. 
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Fig. 5. — Fourier transform of the N-z data for the 
three 6° declination slices on the sky. The black 
solid curve is for the 0° < Dec < 6° slice, the red 
broken curve is for the 40° < Dec < 46° slice and 
the gray solid curve is for the 52° < Dec < 58° 
slice. The numbers are the positions of the peaks 
in the Fourier spectrum when all the SDSS data 
are taken together. See Fig. 7. 
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Fig. 6. — The N-z relation (solid black curve) 
from all the SDSS galaxy data z < 0.35, with red- 
shift bins Sz ~ 10~^. The dot-dashed (red) curve 
is a polynomial fit for the survey function. 
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Fig. 7. — The Fourier spectrum (filled circles, dot- 
ted line) of the SDSS galaxy N-z data (curve in 
Fig. 6) where the DC Fourier component has been 
removed. The sum of 4 Gaussians has been fitted 
to the spectrum (solid black line). 
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Fig. 8. — N-z data (with bins Sz = lO^'^) from 
both regions of the sky. The black curve is from 
the larger region of the sky with 100° < RA < 
270° and gray curve is from the smaller region of 
the sky with 290° < RA < 50°. 



6 but ignores the peaks (the dot-dashed curve in 
Fig. 6) , white noise of the same magnitude as the 
variations of N{z) at higher redshifts was added. 
Its Fourier spectra was then calculated and com- 
pared with the original in Fig. 10. This result 
indicates, as expected, that random white noise 
produces no peaks in the Fourier spectrum, but 
also that the 4 prominent peaks are anomalous in 
relation to what might be expected from a truly 
random distribution of galaxies. It is well known 
that galaxies cluster and also form super clusters, 
long filaments and voids, so some regularity could 
arise from this. From Fig. 10 it is apparent that 
the peaks are significant to at least 4(t, when com- 
pared to the simulated white noise spectrum. 

Table I: Redshift (Az) and distance intervals 
(Ar) for Hq = 100 km,s'^Mpc-^ from N{z) 



Survey 


Peak 


Az 


error 


Ar 


error 


SDSS 


1 


0.0102 


0.0009 


30.6 


2.7 


Fig. 7 


2 


0.0170 


0.0012 


51.0 


3.6 




3 


0.0246 


0.0023 


73.8 


6.9 




4 


0.0448 


0.0131 


134 


39 


2dFGRS 


1 


0.0109 


0.0008 


32.7 


2.4 


Fig. 12 


2 


0.0145 


0.0004 


43.5 


1.2 




3 


0.0244 


0.0035 


73.2 


10.5 




4 


0.0415 


0.0082 


125 


25 



The whole SDSS galaxy N-z data Fourier spec- 
trum is shown in Fig. 7 with the DC component 
removed. The sum of 4 Gaussian curves was fit- 
ted to this. The best fit peak Az values (rep- 
resenting preferred redshift intervals) and their 
standard deviations are shown in Table I. Using 
Hq = 100 kms~^Mpc~^ the preferred redshift in- 
tervals (numbered peaks in Figs 7 and 12) were 
converted to co-moving distances and are listed in 
Table 1. 

2.3. 2D: SDSS and 2dF GRS compared 

In this section the SDSS data are compared 
with those from the 2dF GRS. Reliable redshifts 
were obtained for 229,193 galaxies. An N-z re- 
lation calculated as a function of re dshift z was 
produced in Fig. 17 of CoUess et al hOOlh. This 
is reproduced in Fig. 11 from 2dF GRS data. The 
difference however is that the SDSS data cover a 
larger region of the sky than the narrow 2 degrees 
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Fig. 9. — Fourier transform of the N-z relation 
from Fig. 8. The sohd (black) curve 1 is N-z for 
the whole SDSS data set, the dashed (gray) curve 
2 is for the larger region of the sky with 100° < 
RA < 270° and the dotted (red) curve 3 is for the 
smaller region of the sky with 290° < RA < 50°, 
multiplied by a factor of 3. 
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Fig. 10. — The black curve is the discrete Fourier 
spectrum of the SDSS galaxy N-z data (same as 
the solid curve in Fig. 7) plotted as a function 
Fourier frequency converted to Az interval. The 
gray curve is the Fourier spectrum of the white 
noise enhanced SDSS galaxy N-z data drawn from 
the polynomial fit in Fig. 6. 



of the 2dF GRS data. 

Nevertheless the Fourier spectrum of its N{z) 
data shows very similar peaks to that of the SDSS 
data. See Fig. 12 where the peaks have been fit- 
ted with Gaussians, and Fig. 13 where the Fourier 
spectra from the SDSS and 2dF GRS data are 
compared. The values for the numbered peaks, 
representing the preferred redshift intervals (Az), 
are listed in Table I, along with their standard er- 
rors. Fig. 13 shows the Fourier spectra from the 
whole SDSS data set (solid black curve), the whole 
2dF GRS data set (dotted curve) and a subset of 
49,044 SDSS data where the declination angle is 
within two degrees of the celestial equator (dashed 
curve). The map of the data for the latter closely 
resembles that of the 2dF GRS data as both are 
sampled close to the equatorial plane, and over the 
same right ascension. 

From Fig. 12 it is clear that peak number 2 
is not very significant in 2D samples. Its magni- 
tude is similar to the noise at the higher Fourier 
frequencies. Notice though that the Fourier spec- 
trum of the 2dF GRS data (dotted curve) is similar 
to the narrow equatorial slice from the SDSS data 
(dashed curve). Also the real space map of 2dF 
GRS data seen in Fig. 18 of CoUess et al (|200ll ) 
is very similar to Figs 1 and 2. Both show similar 
arcs and filament structures. 

However peaks 1, 3 and 4 in Fig. 12 are all very 
significant. Also from a comparison of the data in 
Table I, it is clear that the values of the peaks 

I, 3 and 4 are consistent with the respective val- 
ues from the other data set, within their standard 
errors. 

Fig. 13 shows the Fourier spectra from the 
two surveys and except for peak number 2, they 
overlap. The height of peak number 3 from the 
SDSS was normalized to match peak number 3 of 
the 2dF GRS, so they could be compared. It is 
observed that the Gaussian fits are a very good 
match. Peak 3 is at a redshift interval Az = 
0.0246, which is most readily visible by eye in Fig. 

II. Peaks at number 4 have a large height differ- 
ence but similar line shapes. 

Peak number 2 is strongly present in the Fourier 
spectrum from the whole SDSS data set but not 
significant in the 2dF GRS data. The Fourier 
spectrum of the thin equatorial slice from SDSS 
(dashed curve in Fig. 13) is similar to that of 
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Fig. 11. — The N-z relation, with redshift bins 
5z = 10"'^, generated from the data for 229,193 
galaxies in the 2dF G RS data release (see Fig. 17 
of CoUess et al (|2001l) ). The dashed (red) curve is 
a polynomial fit for the survey function. 
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Fig. 12. — Fourier transform of the N-z data from 
the 2dF GRS (Fig. 11) with DC component re- 
moved. The four lowest Fourier frequency peaks 
are shown. The peak data are shown in Table I. 



the 2dF GRS result. Peak number 2 seems to be 
very weak in both samples taken near the celestial 
equator. This indicates that it is not due to insuf- 
ficient statistics. Peak number 1 then represents 
Az = 0.0102, which is also significant in the thin 
equatorial slice. 

3. Mass density fluctuations 

Another way to analyze this data is to calcu- 
late the mass density fluctuations as a function 
of radial distance from the center of the spherical 
distribution. The rms density fluctuations can be 
written in terms of the power spectrum |(5fcp, 



^) ={S{x)5{x)), 



(3) 



whe re (• • ■) indicates and average over all space. 
See (|Kolb fc Turner 1990l) pp. 329-336. The con- 
tribution to {5p/ from a given logarithmic in- 
terval in k is 



A^(fc).l^ 



(4) 



which is the variance or the fluctuations in power 
per logarithmic fc-interval. The rms mass fluctu- 
ations {5M/M) on a given mass scale are related 
to this. 

We use spherical coordinates and calculate the 
average mass density in a spherical volume of ra- 
dius r = ndr where n is an integer and Sr — 
3h~^Mpc. This involves the reasonable assump- 
tion that the number density of galaxies is pro- 
portional to the mass density p{x). Therefore for 
a given radius r 



6M 



(5) 



Below we calculate the averaged density fluctu- 
ations as a function of radius r where 5p is the 
deviation from the average density p, 



AN N{r) 
'Kv' V{r) 



1, 



(6) 



where AiV is the number of galaxies within the 
rtth incremental volume AV determined by taking 
successive radius increments Sr. N(r) is the total 
number of galaxies within the volume to r and 
V{r) is the total volume to r. {6p/p)'^ from ^ is 
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the variance A^(r) evaluated at r. From ([6]) we 
get 



AN AnSrJ27=i{iSry 



Aflr^Sr 



N{r) 



1, 



(7) 



where the subscript has been dropped, Af2 = 
sinOAOAcf) is the sohd angle subtended by the nth 
volume element, which is a constant for data sam- 
pled between fixed RA and Dec angles. Therefore, 



N(n) n2 



1, 



(8) 



where r = n5r has been used, and N(n) is now a 
function of the nth increment. 

Since z — (Ho/cjr it is equally valid to use red- 
shift intervals dz where z = nSz. For the radial 
increment chosen here 5r = 3 h~^Mpc the redshift 
increment Sz — 10^^. It follows from (El) that 



AN{n) (l + n)(l + 2n) 



ELlAiV(^ 



6n 



-1, (9) 



where AN{z) is now the number of galaxies in the 
nth redshift interval, where z = nbz. 

The resulting density fluctuations (bpj p) have 
been plotted in Fig. 14 as a function of redshift for 
both the SDSS and 2dF GRS data. The square of 
these plots give you the variance ^(z) from which 
the power spectrum maybe derived with ap- 
plication of a windowing function and integrating 
over a logarithmic-fc interval. Here we calculate 
P{z) in redshift units instead from 



(10) 



with constants set to unity and where W(z) is the 
top-hat windowing function. 



W{z) 



1, z < zo 

0, Z > Zq. 



(11) 



The Fourier transforms of {&p/p) and [Sp/p)"^ 
were calculated using ([T]) with Us = {6p/p) and 
{5p/p)^, respectively. Fig. 15 shows the results 
for the SDSS data and Fig. 16 for the 2dF GRS 
data. The labeled peaks are listed in Table II, 
where they have been calculated from Gaussian 
fits to the {Sp/p) curves. With the exception of 
Peak 4, the rest are all in good agreement with 



0.006 



0.004 - 



CD 

o 



0.002 




5 10-^ 3 10-^ 10-^ eio-* 4 10-' 

Fourier frequency (Az) 

Fig. 13. — Fourier transforms of the normalized 
N-z data compared. The solid curve is from all 
the SDSS (Fig. 7), the dashed (green) curve from 
a 2° slice through the equatorial plane from SDSS 
and the dotted curve is from the 2dF GRS (Fig. 
11). For the comparison the FTT spectrum from 
the SDSS has been normalized to the peak height 
of peak number 3 in the 2dF GRS result. 
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Fig. 14. — Sp/p plotted as a function of redshift, 
z. The solid (black) curve is for all the SDSS data 
and the broken (red) curve for all the 2dF GRS 
data. 
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Fourier frequency (Az) 



2 10"^ 



Fig. 15. — The solid (black) curve is the squared 
Fourier transform of the SDSS (Sp/p) data from 
Fig. 14, plotted as a function of the redshift in- 
terval (Az). The broken (red) curve is the Fourier 
transform of {5p/ pY derived from the same SDSS 
data. The most prominent peaks are labeled and 
correspond to the labels in Table I and II. 




5 10' 3 10"' lO"' 6 10M10'' 2 10"' 

Fourier frequency (Az) 

Fig. 16. — The solid (black) curve is the squared 
Fourier transform of the 2dF GRS (Sp/p) data 
from Fig. 14, plotted as a function of the red- 
shift interval (Az). The broken (red) curve is 
the Fourier transform of {Sp/ p)'^ derived from the 
same 2dF GRS data. The most prominent peaks 
are labeled and correspond to the labels in Table 
I and II. 



peaks determined from FFT analysis of the N{z) 
counts in the previous section. 

The power spectra P{z) from all the SDSS and 
2dF GRS redshift data are shown in Fig. 17. Both 
can be closely approximated by 11.462 Az. Both 
power spectra deviate from the linear dependence 
and this can be seen in Fig. 18 where this linear 
dependence has been removed. The bottom (red) 
dotted curve is the difference between the two oth- 
ers. There is a large peak near Azo — 0.026 and an 
indication of another near Azq = 0.01, indicated 
by the arrows. 

4. Redshift separation between galaxies 

Another way to determine if galaxies are sepa- 
rated by a periodic spatial interval is to build his- 
tograms by binning the number of pairs of galaxies 
{Npairs) that have the same separation in redshift 
space (Az). Then look for over abundance peaks 
in the histograms. Since redshifts are measured 
radially from the observer at the center of the dis- 
tribution this method detects spatial separation 
with respect to that symmetry. 

Three histograms were built and are shown in 
Fig. 19: 

1. all the 427,513 SDSS data from DR5 release, 

2. the slice of 49,044 SDSS data taken from 
within ±2° of the celestial equator, and, 

3. aU the 229,193 2dFGRS data taken within 
2° of the celestial equator with reliable red- 
shifts. 

In Fig. 19 one can see slight ripples in these 
curves which represent above average correlations. 
Clearly this method is not very sensitive to dis- 
tinguish preferred separations against the back- 
ground of many pairs that have no such correla- 
tion. However one peak is clearly visible as indi- 
cated by the arrow; this is at Az « 0.025. 

By subtracting off a second order (quadratic) 
polynomial from these curves one can separate off 
the background, enhance the weak features and see 
the peaks more clearly, as shown in Figs 20 and 
21. Nevertheless they are still quite weak. In Fig. 
20 the solid (black) curve for Npairs from all the 
SDSS data shows peaks at Az — 0.0103 in good 
agreement with the previous analyses, at Az — 
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Table II: Redshift (Az) and distance intervals 
(Ar) for Ho — 100 kms'^Mpc"^ from (Sp/p) 



Survey 


Peak 


Az 


error 


Ar 


error 


SDSS 


1 


0.0101 


0.0004 


30.3 


1.2 


Fig. 15 


2 


0.0143 


0.0036 


42.9 


10.8 




3 


0.0244 


0.0044 


73.2 


13.2 




4 


0.0557 


0.0130 


167 


39 


2dFGRS 


1 


0.0112 


0.0008 


33.6 


2.4 


Fig. 16 


2 


0.0147 


0.0012 


44.1 


3.6 




3 


0.0251 


0.0020 


75.6 


6.0 




4 


0.0502 


0.0110 


151 


33 




0.06 0.05 0.04 0.03 0.02 0.01 

Fourier frequency (Az) 

Fig. 18. — Power spectra of Fig. 17 where the dot- 
ted (red) linear fit in Fig. 17 has been subtracted 
from both. The solid (blue) curve is for the SDSS 
data and the dot-dashed (black) curve for the 2dF 
GRS data. The dotted (red) curve is the difference 
between these two curves. 




Fig. 17. — Power spectra calculated from SDSS 
and 2dFGRS data. The dotted (red) curve is the 
fitted curve 9.012 x 10"^ + 11.462 Az. 



Fig. 19. — The solid (black) curve (1) is the his- 
togram for all the SDSS data, the dotted (red) 
curve (2) is the histogram for the SDSS slice and 
the dot-dashed (blue) curve (3) is the histogram 
for all the 2dFGRS data. The arrow indicates a 
peak at Az « 0.025. 
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0.0200 and 0.0501 in reasonably agreement, but 
the others are not in agreement. The dotted (red) 
curve from the SDSS equatorial shoe is quite noisy 
without any real clear peaks except the one near 
A2; = 0.011 as indicated. It is possibly the only 
reliable peak considering the subtraction process. 

Fig. 21 shows the N pairs results from all the 
2dF GRS data after a quadratic curve has also 
been removed to enhance the fine detail of the 
peaks. That is shown in the upper curve. The 
lower curve results after a further linear curve was 
subtracted to enhance the first peak. The results 
from the 2dF GRS Npairs analysis gives better re- 
sults than that for the SDSS data. In the former 
case we see peaks around Az = 0.012, 0.025 and 
0.047 in good agreement with the previous analy- 
ses. 



Fig. 20. — The solid (black) curve is the histogram 
for all the SDSS data, the dotted (red) curve is the 
histogram for the SDSS slice after 2nd order poly- 
nomials have been subtracted to show fine detail 
of the peaks. The Az values for the peaks are 
indicated. 
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Fig. 21. — The upper curve is the histogram for 
all the 2dF GRS data after a 2nd order polyno- 
mial has been subtracted to show fine detail of 
the peaks. The inset lower curve is the same after 
an additional linear curve was subtracted to show 
the detail of the first peak. The Az values for the 
peaks are indicated. 



5. Correlation between two surveys 

The correlation between the N{z) counts from 
the two surveys was evaluated after subtracting off 
the background galaxy survey functions, indicated 
by polynomial fits in Figs 6 and 11. This leaves 
only the abundance peaks as a function of red- 
shift which were used in the following correlation 
function. 



{nY.XY~{Y^X){Y^Y)f 



(12) 

where X and Y represent N{z) from the two sur- 
veys and n is the number of AA^^ bins in the N{z) 
distribution-the same number and bin width for 
each survey. 

The function E? is calculated while incremen- 
tally shifting the bin number i, assigned to ATV^ 
in either the X or y distributions, by unity up to 
half the total number of bins. For strong corre- 
lation the unshifted bins (0) value for E? will be 
maximum (where < 1). The result is shown in 
Fig. 22, where the zero shifted value is indicated 
by the arrow. The unshifted value of is max- 
imum except for the the single large peak to its 
right, which can be explained. 

If both distributions are periodic with the same 
period we'd expect also periodic structure in the 
correlation function i?^, which is what we in fact 
see. The Fourier transform of R^ has a single peak 
where Az « 0.027. The strong peak on the right of 
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the central 0-shifted value results when the other 
of the two central abundance peaks in Fig. 11 
aligns with strong peak of Fig. 6. 

6. Discussion 

This analysis finds redshift spacings of Az = 
0.0102, 0.0170, 0.0246, and 0.0448 in the SDSS 
and very similar results from 2dF GRS. The most 
significant peaks in the Fourier spectra, calculated 
from the galaxy N-z data, represent the first, third 
and forth spacings. By combining results from 
both surveys we get the most probable values for 
regular co-moving distance spacings of 31.7 ± 1.8 
h-^Mpc, 73.4±5.8 h~^Mpc and 127±21 h-^Mpc. 
Because the 2dF GRS sampled close to the celes- 
tial equator and the second peak is not prominent 
in the SDSS data when sampled near the celestial 
equator (see dashed curve in Fig. 13), the second 
peak must be related to some redshift symmetry at 
higher declinations. It is strongly seen in the SDSS 
52° < Dec < 58° slice. See peak 2 in the solid gray 
curve m Fig. 5. Peak number 4 (127±21 /i-^Mpc) 
is i n agreemen t with the results of Broadhurst et 



al (Il988l : ll99nh 



The mass density as a function of redshift has 
been calculated from the available data from both 
surveys and the Fourier transforms calculated. 
The results are in good agreement with the Fourier 
transform of the N-z data. Compare Tables I and 
II. The power spectrum is not very sensitive to 
these effects but nevertheless peaks are detected at 
the same redshift periods. Also pair-counts were 
done on both surveys and the Npairs histograms 
showed (albeit weakly) some peaks. The best re- 
sults were gained for the 2dF GRS data. However 
in the SDSS data that method is not sufficiently 
sensitive to detect the effect sought. 

7. Conclusion 

A lot of work has been done in the last 20 years 
trying to show that the spatial structure deter- 
mined from galaxy survey data is in agreement 
with what is predicted by the standard model. 
The standard model is assumed to be correct 
and when evidence question ing t hat conclusion 
is found (Broadhurst et al. ( 1990f ) and evidence 
for concentric shells) a special effort was immedi- 
ately made to show how i t could still be explained 
in the standard model (Kaiser fc Peacock 1991 



Praton et al. 19971 ). But Kaiser and Peacock also 
concluded that their results "did not prove that 
unconventional pictures for large-scale structures 
were ruled out." The appearance of concentric 
rings, like the picket fence, can be explained if ap- 
propriate assumptions are made (e.g. the Bull's 
Eye effect), but this too does not mean that that 
explanation is necessarily correct. Furthermore, 
surely one cannot attribute peaks in the power 
spectrum to a windowing function if the structures 
producing these peaks are already visible in the 
raw data, before the Fourier transform is taken. 

This paper uses a simple Fourier analysis of 
number counts of galaxy redshifts in the SDSS and 
2df GRS galaxy samples to search for periodic red- 
shift spacings. Several spacings are found and are 
also confirmed using mass density fluctuations and 
the redshift separations between galaxies, though 
the latter are less sensitive methods. This analysis 
is similar in many ways to recent analyses carried 
out by other investigators, such as Tegmark et al. 
(j2004l ) where the power spectra obtained from the 
SDSS and 2dF GRS samples are also found to have 
similar power peaks. 

Acoustic oscillations in the plasma during the 
pre-recombination epoch are expected to have 
been imprinted on the distribution of matter in the 
Universe and to have survived until the present. In 
standard model the size, distribution, and cluster- 
ing of galaxies were set in the very early Universe. 
In addition to the windowing function, non-linear 
evolution of structure and scale dependent bias, 
also play a role in our correct understanding of 
this. In this paper we have introduced another 
possible uncertainty into these studies that has 
not previously been considered. We suggest that 
the power peaks we find are also consistent with 
oscillations in the expansion rate at recent epochs. 
As unlikely as this may seem, it is not possible to 
rule it out completely since there is visible evi- 
dence in the raw data for an apparent concentric 
shell structure centered on the observer. It would 
be easier to rule it out if one could argue that 
such oscillations in the Hubble flow are impossi- 
ble. However, astronomers are now arguing for a 
speed-up in the Hubble flow at the present epoch. 
Also a model was recently published (Hirano et al. 
(Hooi)) describing the origin of such oscillations. 

One of the arguments against the power peaks 
being related to oscillations in the Hubble flow 
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is the fact that the P{k) distribution is a rea- 
sonable fit to what is predicted for the standard 
model. However, this does not rule oscillations 
out also. If there were no density fluctuations 
(no clustering etc.), and the Universe was made 
up of a completely uniform distribution of galax- 
ies, oscillations in the raw Hubble flow data would 
likely result in the appearance of concentric shells 
of apparent low and high density surrounding the 
observer because of the light travel time. To de- 
tect this there may be no need to use a windowing 
function. However, even if oscillations are present 
in the low-redshift Universe, any accurate model 
would also have to account for the density fluctu- 
ations (due to clusters, super clusters, etc.) that 
are also clearly visible. 

And some have argued that the peaks found 
in this paper are just the sc al e sizes of clus- 
tering (iKaiser fc Peacock 199ll ; IMo et al. 1991 
Gonzalez et al. 2000f): this was considered im- 
probable by (jYoshida et al 20011 ) but a similar 
analysis needs to be performed on a larger data 
set in order to test this. Nevertheless the meth- 
ods applied to the data here are consistent with 
the Universe undergoing oscillations in its expan- 
sion rate over past epochs. This means in redshift 
space there are preferred redshifts where galaxies 
are more dense due to a slower expansion rate at 
those epochs. 
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Fig. 22. — The correlation function B? between 
the abundance distributions of SDSS and 2dF 
GRS. The zero shifted correlation coefficient is 
shown by the arrow. 
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